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The non-Hermitian quantum dynamics of excitonic energy transfer in photosynthetic systems 
is investigated using a dissipative two-level dimer model. The approach is based on the Green's 
function formalism which permits consideration of decoherence and intersite transfer processes on 
comparable terms. The results indicate a combination of coherent and incoherent behavior at higher 
temperatures with the possibility of exceptional points occurring at the coherent-incoherent crossover 
regime at critical temperatures. When each dimer site is coupled equally to the environmental sources 
of dissipation, the excitonic wavepacket evolves with time with a coherent component, which can be 
^SJ , attributed to the indistinguishability of the sources of dissipation. The time evolution characteristics 

T— I ■ of the B850 Bchls dimer system is analysed using typical parameter estimates in photosynthetic 

■ systems, and the quantum brachistochrone passage times are obtained for a range of parameters. 
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CN ■ I. INTRODUCTION 

■ The theory of excitonic energy transfer has been a topic of interest over several decades in light harvesting systems 
Q-i (LHS) [lHT5[. The Fenna-Matthews-Olson (FMO) complex of the green sulfur bacteria constitutes the 
^l* ■ prototypical LHS for modelling photosynthetic activities [1, [l^, HH, mainly due to its generic features of energy 
Ph ] transfer also seen in other photosynthetic systems [^-01 • In the FMO complex, two types of light-harvesting molecular 
j3 ■ systems complexes, known as LHl and LH2, perform different roles. The LHl is directly linked with the reaction 
qh center (RC) unlike the LH2 complex which interacts with the RC via LHl. The main attraction in light harvesting 

I— I, systems is the exceptionally high efficiencies at which excitation propagates between the light harvesting complexes 
before reaching the reaction center (RC) pigment protein complex Infact, the quantum efficiency in transfer 

, at low illumination intensities reaches close to a value of unity. The importance of mimicking the transfer processes 
^ ■ seen in LHS has obvious applications in artificial light harvesting system such as Ruthenium based complexes [22 24]. 

\ To this end, accurate knowledge of the underlying mechanisms of energy transport and the process by which excess 

. energy in excited states is transferred to their surroundings in LHS is critical to the development of potential sources 
l/^ ■ of energy generation [25l - l28j . 

, Early experimental work in 1960 by Chance and Nishimura [4] showed that some level of photosynthetic activity 

■ occurred in chlorophyll systems close to room temperatures (300 K). The occurrence of quantum coherence during 
photosynthesis was not obvious due to the level of precision of the apparatus employed at that time. However recent 
progress in experimental techniques such as two-dimensional Fourier transform electronic spectroscopy (2DFTES) 
|16,. ,18;, .29, ,30j has confirmed that quantum coherence is indeed conserved for relatively long times (up to 1 picosecond) 
in excited states. Using 2DFTES [Tgl, the FMO complex of green sulfur bacteria showed a surprisingly long coherence 
time of about 700 fs at 77 K, and also at room temperatures (277 K) in a related work [3l|. These ambient temperatures 
are generally considered adverse for quantum features such as superpositions and non-classical correlations to be 
maintained for a reasonable period of time. However the experimental results [l^, [l^, [29l - l3]| appear to suggest the 
critical role played by non-local quantum effects in energy transfer mechanisms in LHS. To this end, the application of 
advanced measurement tools involving spectrally resolved, 4-wave mixing measurements [32| is expected to introduce 
greater depth to the study of quantum correlations in photosynthetic systems which undergo both decoherence and 
dissipation due to interactions with the environment. 

Calve et. al. (ssj showed recently that for a system of two interacting and parametrically driven harmonic oscil- 
lators undergoing dissipation, entanglement can exists at any temperatures. Although this specific quantum system 
appear distinct from quantum models which describe LHS, the general consensus is that the slower processing speed 
inherent in classical models is inadequate for LHS compared to the processing attributes of quantum systems. This 
viewpoint has been the focus of many recent works [9l. [3#j40| involving LHS. A variety of methods involving the Red- 
field approach employing the Born-Markov and secular approximations [4llj43t. the non-Markovian approach which 
involves solving the integrodifferential equation using perturbation theory [44], and recently the more sophisticated 
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and numerically intensive reduced hierarchy equation approach [455 have been employed to examine the propagation 
of excitation in LHS. In light-harvesting complexes with strong system-bath coupling, a generalized BlochRedfield 
(GBR) equation approach was used to show that optimal efficiency in energy transfer cannot be optimized with 
respect to the temperature and spatialtemporal correlations in noise [46^. Thus classical approaches to the excitonic 
transport which involve the hop ping model [2, Q have given way to more sophisticated quantum mechanical models 
of exciton propagation jol. [ill I45l - l47l | which incorporate quantum coherence aspects during the energy transfer process. 

Greens functions are known to provide an effective description of the quantum evolution of non equilibrium quantum 
processes [49-53], and here we use this approach to examine the influence of dissipative processes on the dynamical 
evolution of the exciton population in LHS. While several works have considered coherence and bath correlations 
[3^ [35I I4H as possible factors for the high efficiencies of LHS, here we consider two features which may account for 
the quantum properties in LHS: (a) the persistence of coherence when the dissipative coupling of each monomer to its 
surrounding environment becomes equivalent to that of the adjacent monomer, and (b) the appearance of exceptional 
points at critical temperatures. The simple model of the dimer is used to show that in the event that the coupling of 
each subsystem to the phonon bath or other sources of dissipation is equal, the overall time evolution still retains a 
coherent component, an idea that was first proposed by Stafford and Barrett [65| in the context of the decay of super 
deformed nuclei systems. This key feature implies that a degree of coherence is maintained even at temperatures that 
would otherwise be considered adverse for any superpositions of the quantum states. Accordingly we consider the 
maintenance of coherence, as quantified by the time duration of oscillations in the population difference at the two 
site model, as a measure of the photosynthetic efficiency. In general, the photon is absorbed by a network system of 
pigments and transported to the reaction center (RC) trap so that efficiency is based on the success of being trapped at 
the RC. The connection between quantum coherence and photosynthetic efficiency is justified as the loss of coherence 
disrupts the likelihood of the propagating exciton of reaching its destination due to decoherence processes. 

Unlike earlier works on LHS, we take advantage of the non-Hermitian features inherent in any open quantum system 
which encounter dissipative forces. We also aim to examine the underlying factors which are linked with the non- 
Hermitian features and which may account for the observed long-lived coherence in LHS. Open quantum systems with 
non-Hermitian components evolve in ways which are vastly different from quantum systems associated with a purely 
Hermitian Hamiltonian. For instance, the states associated with an Hermitian Hamiltonian have long lifetimes, while 
those of the non-Hermitian Hamiltonian have a finite lifetime. Non-Hermitian resonances are associated with complex 
eigenvalues with real (imaginary) component that yield the energy (resonance width), and possess non-orthogonal 
eigenvectors. Thus states which are orthogonal under the ordinary inner product in the Hermitian quantum space 
are allowed to be non-orthogonal in the non-Hermitian case. The real state energies of a Hermitian Hamiltonian give 
rise to the avoided level crossing [s^ l when a single parameter is varied, which is not necessarily true in the case of 
a non-Hermitian Hamiltonian. The presence of non-Hermitian terms is critical to the occurrence of dynamical phase 
transitions as well. Moreover the appearance of degeneracies such as exceptional points (EPs) is a unique feature of 
non-Hermitian quantum systems, and may assist in distinguishing classical and quantum modes of transport in LHS. 

A topological defect like the exceptional point occurs when two eigenvalues of an operator coalesce as a result 
of change in selected system parameters, and the two mutually orthogonal states merge into one self-orthogonal 
state, resulting in a singularity in the spectrum ^55]. The critical parameter values at which the singularity appear 
are called exceptional points. A notable feature associated with exceptional points is the violation of the adiabatic 
theorem, that is the switching of an initial eigenstate to another eigenstate when the system parameters are altered 
adiabatically around the exceptional point. In contrast, Hermitian systems exhibit only berry phases ^56i], with 
no change in the initial eigenstate as the system parameters are altered adiabatically. Thus no exceptional points 
can exist in Hermitian systems and at degenerate points the eigenstates exist only in a two-dimensional subspace 
spanned^by vectors. EPs have been observed in experiments involving microwave billiards [57l - [59| . semiconductors 
cavities 60] and photo-dissociated vibronic resonance states of the H^ molecular ion (6]| . Exceptional points may be 
associated with the transition or crossover points at which coherent to incoherent tunnelling occur in open quantum 
systems. The detection of the exceptional points are still under active investigation [g^. One possible technique 
involves the detection of the switching process, in which two states around an exceptional point swap states and thus 
behave distinctly from other states which remain unchanged when the selected parameter is changed adiabatically 
[ssj . However the experimental detection of such topological changes still remains a challenge. 

This work is organized as follows: In Sec. (|TT]) we present the theory of the exciton transfer using Green's function 
formalism and analyse the infiuence of the intersite coupling energy and dissipation rates on the coherence properties 
of the exciton in a simple dimer system. In Sec. (IIIII) . we examine the photosynthetic qubit system and discuss 
the conditions under which exceptional points appear. In Sec. (jIVp . the dimer model is applied to the B850 Bchls 
system and numerical estimates of the time scales of coherent oscillations in the exciton population difference are 
obtained. Estimates of the critical temperatures at which exceptional points occur are also evaluated for a range of 
environmental dissipation strength differences. Lastly, in Sec. (jVj, the quantum brachistochrone passage times are 
obtained for a range of parameters in photosynthetic systems. A brief discussion and conclusion are also provided in 
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Sec. m 



II. EXCITON TRANSFER USING GREEN'S FUNCTION FORMALISM 

The Hamiltonian describing the propagation of the exciton in molecular systems is based on the tight-binding model 



Hex — 



(1) 



where Bj (Bj) is the creator (annihilation) exciton operator at site j. AE, the on-site excitation energy and Di,m is 
the dispersive interaction matrix element which determines the energy difference between a pair of excited electron 
and hole at a molecular site and ground state electrons at neighboring sites ^2] . V the exciton transfer matrix element 
between molecular sites at / and m. In the case of the dimer system with just two coupled sites (labeled / and m), 
Eq. ([T]) is greatly simplified 



Hex — EfBjBi + E^j^Bjj^B,, 



vb\b„ 



(2) 



where Ei,Ejn are the exciton energies at sites I, m and the subscript denotes the absence of lattice site fluctuations. 

In the presence of a non-Hermitian decay terms and phonon bath reservoirs at each site, the exciton dynamics of 
the dimer system is determined by the following Hamiltonian 



Hd 
Hp 

Hpri 
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where the isolated exciton Hamiltonian, Hex is given in Eq. ([2]). We allow the intersite tunnelling amplitude V to be 
influenced by a phonon correlated environment via a Franck-Cordon (FC) factor [g^l 



Vr ~V exp 



J(w) 



coth 



duj 



(7) 
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where the spectral density function is given by J{uj)— —J2q l^a(9)P S{uj — ujq). We assume that the correlated exciton- 
phonon interaction term Fa{q) to be the average {^[Fi{q) + ^^(all of the interactions at sites, l,m. The argument 
of the exponential term in Eq. ([7]) is known as the FC factor [6J| and its dependency on the temperature varies 
according to the form chosen for J{uj). 

The non-Hermitian dissipative Hamiltonian in Eq. Hd is associated with the decay rate of the exciton 
at each site j — l,m. For simplicity in numerical analysis, we have incorporated several processes such as exciton 
annihilation due to recombination and trapping effects into Q , with the requirement that the decay rate is dependent 
on the site j. Due to differences in environmental conditions at the two sites of the dimer Q ^ (^m, however the 
possibility that Q « Cm cannot be excluded. Hp denotes the sum of phonon energies at the two sites and b^^ j{bqj) 
is the creation (annihilation) phonon operator associated with the phonon bath at site j, with frequency and 
wavevector q. Hep represents the sum of the exciton-phonon interactions at the two sites, and Fj{q) quantifies the 
interaction strength which is assumed to be linear in the phonon operators and which is dependent only on the phonon 
wavevector q. We consider that both dissipation and dephasing processes occur as a consequence of exciton-phonon 
interactions. The dissipative mechanisms associated with the exciton-phonon interactions introduce a non-Hermitian 
term in the Green's function (see Eq. PH)) below) which account for the irreversible loss in the exciton population. 
While the dephasing term leaves the overall exciton population unchanged, it contributes to incoherence by changing 
the non-diagonal components of the density matrix. 

The Green's function for an exciton at time t which yields the dynamical measure of the amplitude of the electron 
(hole) to propagate forward (backward) in time is given by f53] 



Gl^rait) ^ -m){{Bl{t)Bl}\) 
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where Q(t) denotes the step function and the exciton Green's function is averaged over the quantum mechanical 
properties of the propagating excitonic wavefunction that is dependent on the phonon bath and contribution from 
the dissipative terms. Following the approach by Stafford and Barrett [HHl, the Fourier transform Gi,m{E) is given in 
terms of the energy E as Gi^m{E) = dtGi^m{t) e*^*. For the simple case of the dimer system with two coupled 
sites of energies Ei,Em, the inverse of the Green's function Gi^m{E) is obtained as 



E-E1-A1+ iri+ -Vr 

— Vr E — Em — 



1 

2^- 



(9) 



where ry is a very small number and Aj is the shift in the exciton energy due to its interaction with phonons. This 
polaronic energy shift can be evaluated using standard techniques based on the polaron model [i^. Without loss in 
generality, the polaronic shift term can be incorporated as a constituent of the site energy Ei or Em, and therefore 
we drop this term from now on. E is the self-energy due to dissipative interactions with the phonon bath and other 
processes associated with exciton recombination, annihilation and trapping which is expressed as 



-i{lph,i + ld,i) 






-i{lph,m + ld,m) 



(10) 



Here we distinguish the phonon-related dissipative term ^ph,j from the dissipative term arising from non-phonon 
related recombination and trapping effects which is denoted by ^d.j- lph,j at site j is obtained as 



lph.„ 



N- 



F,{q)\H{hL.- 



E - Ap) 



(11) 



where Em is the mean energy of the exciton band, E is the exciton energy and Ap can be interpreted as the dissipative 
or trap depth [G^j which results in the irreversible loss of the exciton. Both Em and Ap appear as phenomenological 
constants and their associated values are unavailable for the FMO complex. For the purpose of obtaining numerical 
values in the next section, we select a range of values for the cumulative term (jph.j+Jd.j) and examine its influence 
on the coherence properties of the excitonic dimer. 

Using the inversion procedure adopted by Stafford and Barrett [65], we obtain 



G{E) = 



E-Ei + -ji 



E — E„ 



Vr 



E - Em + -^Jr, 
Vr 

E-Ei + f 7i 



Vr 



(12) 



where "fm = lph,m + ld,m and 7/ = ^ph,i + ld,i- We consider that the exciton in the simple dimer system is initially 
localized at site j = I at time t = 0, thus the probability Pu that the exciton remains at site j = I is determined using 
Pii{t) — |Gii(i)|^. Explicit expressions for Pu{t) and the probability that the exciton has propagated to site j = m, 
Pim{t) can be obtained for equal site energies, Ei = Em in the coherent tunneling regime {2Vr > 7*) 



Pu 

Pirn 



n 7^ n 



-7t 



4K- 



(13) 



where fl = (41^^ — 7*^)^^^, 7* = ^(7m — 7i) and 7 = ^(7m + 7;)- While the contribution of the dissipative term, 7 
appears in the exponential function (see Eq. ([l3l) ). the dimer can be seen to still exhibit coherence with an effective 
Rabi oscillation that is influenced by the environment. The system however undergoes incoherent tunneling with loss 
of Rabi oscillations at 2Vr < 7* and we obtain 



Pu = 



Ph 



- .-7* 



cosh —t smh — t 

2 n 2 



\v 



'It 



■ ,2^ 

smh — 



(14) 



Due to the presence of an non-Hermitian term in Eqs.®, (ITU|) . the total probabilities, Pu+Pim < 1 is not conserved 
and there is loss of normalization which is dependent on the dissipative terms, 7; and 7m- 
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We note that in the case of unequal energies, Ei ^ Em-, the dimer system exhibit both incoherent and coherent 
characteristics with a combination of the terms which appear in Eas. (|13p and (|14p. In the event of equivalent dissipative 
couplings at each subsystem, we obtain 7; = 7^, 7*=0, and the imaginary term present in vanishes as noted earlier 
in Ref.[68| for a study involving quantum dot systems. Despite the presence of non-zero 7; and 7™, Rabi oscillation 

of frequency f2 = ^ 414-^ + (Ei — Em)^ occurs. The subsystems of the dimer can be considered to be unmeasured 
by the environmental sources when there is equivalent couplings to the dissipation channels, an observation that was 
first noted by Stafford and Barrett 65:] for super deformed nuclei systems. In this regard the indistinguishability of 
the source of decoherence preserves the coherence of the excitonic dimer in photosynthetic systems. 



W (b) (c) 




FIG. 1: (a) The exciton population difference, AP=Pii — Pi,„ as a function of time t, and the bare intersite coupling energy, V 
at dissipation rates 7m=7i=0.1 for the degenerate case {Ei — Em)- The units are chosen such that h—1, and the phonon bath 
response time, ujq — I. 

(b) AP=Pii — Pim as a function of time t, the bare intersite coupling energy, V at dissipation rates 7,„=0.1, 7!=0.5. 

(c) AP—Pii — Pim as function of time t and dissipation rate 7,„, at 7; =0.2 and V=3- 

In order to examine the effect of the intersite coupling energy, V and the dissipation rates , 7/ on the exciton 
dynamics, we consider the bare intersite energy V without the Franck-Cordon (FC) factor [gJ]. Figure [T^,b,c shows 
the coherence properties as reflected in the exciton population difference, AP—Pu — Pi„i as a function of time and 
the bare intersite coupling energy, V - We note that an increase in the strength of intersite energy V increases the 
time period over which the population difference, AP which is a signature of coherence, is maintained. As illustrated 
in Figure [ij:, the gradual increase in one of the dissipation rate 7™ erodes the coherence in the dimer system as time 
progresses. 



III. PHOTOSYNTHETIC QUBITS AND APPEARANCE OF EXCEPTIONAL POINTS 



Following Ea. p2)l and the form of the occupation probabilities in Eqs. ([T3)) and p4)) . the symmetric and anti- 
symmetric states of the dimer system at the resonance point {Ei=Em) is obtained as 

\Xsit)) = e-^*/2 (^cos^t-icos0sin^t^ |1) 

-l-ie"^*/^ sin sin |m) 

IXaW) = e-^*/2 (^cos^t-icos0sin^tj |1) 



-ie ^*/^sin6'sin^t |m), 



(15) 



where cos0=-^ and VL and 7* are defined below Eqs. ([T3|) . The excitonic qubits states in Eqs. (fT5|) are coded using 
the relative position of the exciton via the basis set, (|1) , |m)) 

|1) = |H,)®|0)„ 
|m) = |0),®|H„,), (16) 

where |Hi) (|Hm)) denote the excitonic state at site I (m), and the state |0); {\Q)m) denote the ground states at site 
I (to) (the absence of exciton). Other than the qubit states, \xs{t)) and \xa.{t))i entangled states of the form |0)( |0)j 
and I Hi) |Hm) may also result in the dimer system. 
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At unequal energies, Ei ^ E„i, the states equivalent to the symmetric and asymmetric states in Eq.([T5]) possess 
eigenenergies of the form 



E±=Ei+ E„, ± {E„, - Elf + AVr\ (17) 

where the complex energy levels Ei—Ei — 17; and Ej^—Em, — i^m and the dissipative terms 7^ and 7^ are specified 
below Eq. p^ . The excitonic qubit oscillates coherently between the two dots with the complex Rabi frequency 

E+ — E-=Cd= 2 [[{Em — Ei)]^ + 4T4^] . By setting Ei = Em, the eigenenergies corresponding to Eq. pSl) can be 
easily evaluated. The real component of uj correspond to Rabi oscillations while the imaginary component correspond 
to the rate of incoherent tunneling [s^l ■ The role of a complex Rabi frequency has been discussed [68] in the context 
of the double quantum dot system coupled to a continuum of states. 

Due to the presence of non-Hermitian terms, we write the adjoint symmetric and anti-symmetric states correspond- 
ing to the states in Ea.([T5l) as 

\Xs{t)) = e-^*/2 (^cos^t + icos0sin^<^ |1) 

-ie~^*/^ sin sin — t |m) 
2 ' ' 

IXaW) = e-^*/2 ^cos^t + icosgsin^tj |1) 



+ie sinSsin |m) , 



(18) 



In the case of strong dissipative processes and at rising temperatures, it is likely that the condition 2Vr — 7* 
will be satisfied and 51=0. This signifies the appearance of the exceptional point (ssj when both coherent and 
incoherent tunneling regimes merge, and the two eigenvalues coalesce to represent just one eigenfunction. We obtain 

Pu = [l- ^) e-^\ Pim = e-^^ and the population difference Pu ~ Pim^ (1 - 7*t) e'^^ 



IV. APPLICATION TO THE B850 BCHLS DIMER MODEL 



The LH2 B850 pigment-protein complex which interacts with the reaction center via LHl, consists of a cyclic array 
of polypeptide heterodimers, with each polypeptide pair of a/3 subunit, containing three Bchls. A pair of Bchls is 
located near the outer membrane surface of the complex, while a single Bchl is present near the inner membrane as 
detailed in Ref . f48| . A similar cyclic structure (B800 ring), but one with a smaller number of Bchls accounts for the 
800 nm absorption spectrum of the LH2 complex. We utilize a simple model which consists of two coupled B850 
Bchls, each BChl with a two-state structure forming the basis of the Qj^ one-exciton states. The Qj, bandwidth is 
associated with the lowest excited state of the BChl molecule. The two-state model of the BChl interacts with its 
surrounding phonon bath with properties similar to that of a spin-boson system [64| . and forms an integral part of 
light-harvesting energy transfer mechanisms in the LH2 B850 pigment-protein complex. 

We employ the Hamiltonian of the form given in Eq. ([3]) to represent the B850 Bchls dimer, with the intersite 
coupling energy V—2bQ cm~^ and consider the two site energies to be equal, that is Ei « Em- We use the spectral 
density function associated with the overdamped Brownian oscillator (69i] 

J{u)) = -AfcWo — 2 (19) 

where Ab is the bath reorganization energy and wq is a frequency cutoff for the bath. We use values of Xb— 200 
cm^^ and hw^j — 50 cm~^ which are acceptable estimates in photosynthetic systems j^, [l^ [l^, [l^, . It is expected 
that the bath reorganization energy Ab has an implicit dependence on temperature, however here we assume that it 
remains independent of temperature for the sake of obtaining some numerical estimates of the dynamics of energy 
exchange in the B850 Bchls system. 

The effect of temperature on the time evolution of the population difference, AP=P/; — Pimi at various dissipation 
levels is shown in Fig. [5^,b,c, While oscillations are preserved at small dissipation levels up to moderate temperatures 
(80 K), there is less exchange of energy between the two sites at the comparatively larger temperatures (140 K). 
Moreover the probabilities of site occupation Pu and Pim becomes zero due to outflow of energy to the dissipation 
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(a) (b) (c) 




t(ps) t(ps) t(pg) 

FIG. 2: Fading out of the exciton population difference, AP=Pii — Pim, with increase in time (ps) and temperature (in K) at 
various dissipation strengths, 7m=7; = (a) 2.5 cm~^ (b) 5 cm~^ ("ilO cm~^. The bath reorganization energy Ai,=200 cm~^, 
ftwo = 50 cm^^ and bare intersite coupUng energy, l/=250 cm~^ 0, E) E, 13 • The colour shading range from red for the 
maximum population difference of 1 to blue for the minimum population difference of -1. 



(a) (b) fc) 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



t (ps) t (ps) t (ps) 

FIG. 3: Population difference, AP=P;; — Pim as function of time (ps) and temperature (in K) at various dissipation strengths, 
7m=7i = (a) 2.5 cm~^ (b) 5 cm~^ (c) 10 cm~^. The bath reorganization energy is set lower at Ai,=150 cm"'^. All other 
parameters {T>W(i,V) and colour codes are set at the same values as specified in the caption of Fig. [21 



sources. We note increased oscillations in AP as illustrated in Fig. ISt^ibjC when the bath reorganization energy 
is reduced to a lower value of Ab=150 cm^^. As mentioned earlier, the results in these figures are based on the 
assumption that the reorganization energy is independent of the temperature, which has been used in earlier works 
[l^|48|. Currently, the explicit dependence of \b on the temperature is lacking in the literature. The effect of including 
the temperature dependence of Ah is expected to alter the quantitative estimates of the timescale of oscillations by 
a small factor, however we expect the salient qualitative features to be preserved. It is to be noted that the time 
estimates obtained using our simple dimer model is in partial agreement with experimental results [29l [sij which 
show that coherent oscillations persist up to 0.3 ps even at physiological temperatures for a similar photosynthetic 
system. The oscillations in population difference AP as illustrated in Figs.^J |3]are clearly dependent on the strength 
of exciton coupling to the sources of dissipation and the estimates in the range (2 to 10) cm~^ employed here may 
be lower than those in photosynthetic systems. Accordingly, recombination and trapping effects which are common 
sources of dissipation, are expected to have a significant effect on the photosynthetic activity of light harvesting 
systems in general. 

Fig. [3] illustrates the dependence of the critical temperatures (K) which occur at exceptional points, as a function of 
the dissipation strength difference, A7=|7„i — 7;| based on Eg. dTTl) . The results indicate that some degree of coherence 
can be obtained at physiologically higher temperatures (> 70K) as long as the dissipation strength difference, A7 
is kept low. It is to be noted that the estimates of the critical temperatures in Fig. 2] is based on a simple dimer 
model which circumvents comprehensive details that incorporate the energy bias and solvent relaxation rate, various 
kinetic behaviors of the electron transfer mechanisms and dynamic bath effects as considered in Ref.fT^]. In general, 
the Redfield model employed here works well as a weak coupling model, however its validity is very much dependent 
on the partitioning between the exciton Hamiltonian and the intersite coupling and cxciton-phonon terms. In Eq. 
([3]) we have chosen the exciton states as the zeroth order term, which may not provide a reliable description of 
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Ay (cm ^) 



FIG. 4: Critical temperatures (K) at which exceptional points occur, as a function of the dissipation strength difference, 
A7=|7m-7!|. All other parameters {huo, Xb, V) are set at the same values as specified in the caption of Fig. [2] 



the dynamics in light- harvesting complexes, under certain environmental conditions [46|. It has been noted earlier 
that the secular Redfield equation [4^ leads to an overestimation of the environment-assisted transfer process which 
results in inaccuracies in the prediction of the coherent-incoherent transition region. One approach of overcoming this 
problem is by extending Eq. ([3]) to dressed excitons or composite exciton-phonon states with hybrid exciton-phonon 
characteristics. This can be done using suitable canonical transformations of exciton and phonon operators to a 
transformed Hilbert space, however analytical solutions can only be obtained under strict limits. While the analytical 
model employed here offers advantage in the computational speed of prediction of the exception point, we expect some 
quantitative changes (less than a order of magnitude) in the critical temperatures shown in Fig. |4] when improved 
models incorporating more realistic features are utilized. Nevertheless, the qualitative trends, especially changes in the 
exceptional point temperatures with dissipation strength differences, are expected to be unaffected by the assumptions 
inherent in the Redfield model for the range of A7 and bath reorganization energies considered in this work. In this 
regard, the underlying principles associated with the appearances of exceptional points can be generalized to studies 
which include more elaborate system parameters. Hence there is ample scope to extend the intermediate complexity 
of the current work to scrutinize the factors which influence the crossover point from coherent to incoherent dynamics, 
and influence the appearance of exceptional points in more realistic systems with strong system-bath couplings. 

The connection between the exactly solvable two-state model to light-harvesting energy transfer processes should 
also be viewed in the context of the large molecular structure of realistic photosynthetic systems. Extrinsic factors 
such as network size and topological connectivity present in large molecular structures assume vital roles in the 
optimum energy transfer processes, and consequently determine the coherence properties of the entangled exciton 
in large molecular structures. In a recent work [Tlj . the mapping of quantum energy transfer processes to network 
kinetics with nonlocal connectivity was examined in three-level systems, linear-chains, and other closed networks. 
The results indicated that environments can be optimized to maximize energy transfer efficiency, with prediction of 
phase-sensitive interferences in closed-loop configurations. These factors provide sufficient framework to examine the 
appearances of exceptional points within the structure of large topologically connected structures and their possible 
influence on phase-sensitive interferences [7l| . 



V. THE QUANTUM BRACHISTOCHRONE PASSAGE TIMES IN PHOTOSYNTHETIC SYSTEMS. 

The exact role of exceptional points in photosynthetic systems remains unclear, and it is likely that these points 
are linked to the quantum non-Hermitian brachistochrone problem [t^ . pTSj . In the general brachistochrone problem, 
the minimum time taken to transverse the path between two locations of a particle is determined, this becomes the 
quantum brachistochrone problem when the time of evolution of the quantum system between two states is required. 
It has been shown that the passage time of evolution of an initial state into the final state can be made arbitrarily small 
for a time-evolution operator that is non-Hermitian but j-T-svmmetr ic |73ll. This result has recently been generalized 
to non-PT-symmetric dissipative systems [t^. The results in Refs. |73l. |74|| appear to suggest that the timescales of 
propagation in non-Hermitian quantum mechanics are faster than those of Hermitian systems. Using the approach 
detailed in Ref. [T^I, we compute the passage time, Tp taken for the |1} state to make a transition to the |m) state as 
a function of temperature and the dissipation strength difference, A7=|7,„-7;| (Fig. [5]). 

As expected, Tp increases with temperature, however it decreases with the dissipation strength difference, A7. The 
passage times (in ps) obtained here are small compared to typical estimates (f ns) of the fiuorescence decay times of 
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FIG. 5: Passage time from the |1) to |m) states as a function of the dissipation strength difference, A7=|7m-7i| and temperature. 
The bath reorganization energy Ab=150 cm~^, hujo = 50 cm~^ and bare intersite coupling energy, ^=250 cm~^ 

the light harvesting complex [isj . This suggest considerable delocalization, with implications for non-trivial quantum 
correlations amongst various pigment sites. However the explicit link between the entanglement properties and the 
brachistochrone passage times is not well understood. There is some consistency of the results in Fig. [5] with earlier 
results which show loss of coherence at physiologically higher temperatures (Fig. 2]), however the decrease of Tp with 
the dissipation strength difference, A'y=\jrn-li\ may be associated with a decrease in efficiency of the energy transfer 
mechanism. Clearly there are open questions with regard to the entanglement properties of quantum systems which 
obey non-Hermitian quantum dynamics, these properties are comparatively well established in Hermitian systems. 
For instance there is lack of consensus on the properties of the non-Hermitian eigenstates, as while it is known 
that real eigenstates are associated with localized states, eigenvalues of non-Hermitian eigenstates may have complex 
delocalization features which are not well understood i75|. 

The results obtained here can be compared to earlier works (tgI - ItsI I of coherent exciton dynamics which show that 
the environment not only allow the preservation of some coherent features but also assist in the quantum transport 
in light-harvesting systems, even at ambient temperatures. However the role of exceptional points have not been 
implicated in the earlier works, and in this respect, the estimates of the quantum brachistochrone passage times here 
is expected to bring greater depth to existing studies in this field. Hence the solvable model examined in this work 
suggest features that have not been predicted in earlier works. Nevertheless, the role of new measures of quantum 
correlations which possess more generalized properties such as the quantum discord [tqI Isoj in the photosynthetic 
process needs further investigation, including examination of the intricate links between the exceptional point and 
quantum correlation measures. 

Lastly, we have presented calculations associated with a dissipative two-level dimer model which incorporates decay 
terms due to coupling to a macroscopic environment. The results when applied to the B850 Bchls dimer system can 
account partly for the observed high efficiencies of energy transfer in light harvesting systems, and highlight the critical 
link between photosynthetic coherences and dissipation strength differences at the two sites of the dimer. There is 
greater efficiency in energy transfer as long as the distinguishability of the subunits of the dimer is minimized, an 
analogous observation was made for the system of quantum dot systems connected to conducting leads in an earlier 
work'68] . Thus in the presence of a small dissipation strength difference, A7 the excitonic coherences appear insensitive 
to the environmental changes, provided the subunits are not distinguished by varying the strength of dissipation. The 
most important result of this study is however related to the appearance of the exceptional point which has not 
received attention in earlier works on photosynthetic systems. Despite the frequent mention of the crossover point 
from coherent to incoherent dynamics at increasing temperatures in the literature, the presence of the exceptional 
points has been overlooked. The evaluation of the quantum brachistochrone passage times is also unique to this work. 
However the role of exceptional points need further investigations, as there are many unanswered questions as to how 
such points actually related to the high efficiencies measure in light harvesting systems. 

It is to be noted that the model used in this work can be applied to obtain numerical estimates of quantum coherence 
in other photosynthetic systems which possess more realistic spectral functions associated with the environmental 
couplings . The model used here can also be extended to examine the properties of the exceptional points and 
the quantum brachistochrone passage times in a dimer model where the onsite and coupling energies differ as well. 
It is expected that the qualitative features involving the appearances of exceptional points, and shown within the 
intermediate complexity of the current work, are expected to be retained in more sophisticated models. In this regard, 
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the results obtained in the current work provide important guideUnes for future investigations in more sophisticated 
models, and in the detection and control of quantum effects in artificial photosynthetic systems. 
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